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Description 

This invention relates to image processing. It finds particular application in conjunction with image en- 
hancement, image smoothing, and other image improvement techniques for magnetic resonance images and 

5 will be described with particular reference thereto, it is to be appreciated, however, that the present invention 
is also applicable to improving digital x-ray images, computed tomographic images, nuclear camera images, 
positron emission scanners, and the like. 

Techniques capable of producing images of the internal structure of a body, such as are required for ex- 
ample, for medical diagnostic purposes, have commonly been subject to image degradation from noise, system 

10 imperfections, and the like. Various image processing techniques have been utilized to remove the effects of 
the noise. See for example, "Digital Image Enhancement A survey" Wang, et al., Computer vision, Graphics, 
and Image Processing, Vol. 24, pages 363-381 (1983). In one technique, each pixel was adjusted in accordance 
with the mean of surrounding pixels and the variance or difference between the surrounding pixels. Each filter 
enhanced pixel value g'(i j) was a weighted average of the local mean and variance values: 

« _ 9'P = 9(iJ) + k[g(ij)-g(i,j)] (1), 

where g(i,j) was the local mean, g(i,j)-g(i j) was the variance, and k was a constant that weighted the relative 
contributions therebetween. It is to be appreciated that when k was set larger than 1 , the variance or difference 
between the local mean value, hence the contribution of the measured gray scale level of the pixel (ij) was 
magnified. As k was set smaller, the image was smoothed or blurred as if acted upon by a low-pass filter. At 

20 the extreme at which k was set equal to zero, each pixel value was replaced by the local mean of the neig h boring 
pixel values. 

One of the drawbacks in this technique resided in selecting an appropriate value for the weighting factor 
k. The smaller k was set, the more the image was blurred and the more difficult it became to withdraw accurate 
diagnostic information. As k was set larger, edges and fine details, including noise, became enhanced. Fre- 
25 quently, in a medical image, the selected weighting factor k was too large for some regions and too small for 
other regions. 

"Digital Image Processing by Use of Local Statistics" by J.S. Lee, Naval Research Laboratory, Washington, 
D.C. (1980), recognized that a different weighting factor k could be selected for each pixel to be enhanced. 
Specifically, Lee suggested setting the k for each pixel equal to the square root of the ratio of a preselected 
30 desirable local variance to the actual local variance of the selected pixel. Although the Lee pixel variable weight- 
ing factor achieved better resultant images than the constant weighting factor, there was still room for improve- 
ment. 

It is an object of the present invention to provide a method of and apparatus for imaging utilising a novel 
technique for image improvement. 

35 In accordance with one aspect of the present invention there is provided a method of imaging comprising: 
collecting image data; converting the collected data into an electronic image representation which Includes 
an array of pixel values; for each pixel value determining (i) a variance between the pixel value and neighboring 
pixel values and (ii) an average of the neighboring pixel values; and replacing each pixel value of the array by 
a weighted combination of the said pixel value and the average of pixel values neighboring said pixel value; 

40 characterised in that the method further includes determining noise statistics of the collected data; and in that 
the weighting is in accordance with the determined noise statistics and the variance between said pixel value 
and its neighboring pixel values. 

In accordance with another aspect of the present invention there is provided an imaging apparatus com- 
prising: data generating apparatus for generating image data indicative of at least a selected region of a body; 

45 a data memory means for storing image data from the data generating apparatus; a transform means for con- 
verting the data from the data memory means into an electronic image representation which includes an array 
of pixel values; an average pixel value determining means for determining an average pixel value of pixel values 
neighboring each pixel value; and combining means for adjusting each pixel value, in accordance with the cor- 
responding said average pixel value; characterised in that the apparatus further includes a data noise statistic 

go determining means for determining a noise statistic from the data stored in said data memory means; and in 
that said combining means combines each pixel value with the corresponding average pixel value and the data 
noise statistic. 

A first advantage of the present invention is that it provides a highest appropriate filtering level for each 
individual image. Because the filter function is image noise dependent, more noise free images are not over- 
55 filtered or blurred unnecessarily. 

Another advantage of the present invention is that it is regionally adaptive. That is, noisier regions of the 
image are filtered more than regions of each image with less noise. 

Other advantages of the present invention reside in an improved technique for determining image wide 
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noise, in achieving an improved signal-to-noise ratio without a loss of spatial resolution, and in an improved 
technique for determining image noise before image reconstruction. 

Still further advantages of the present invention will become apparent to those of ordinary skill in the art 
upon reading and understanding the following detailed description of the preferred embodiment. 

5 One method and apparatus in accordance with the invention will now be described, by way of example, 

with reference to the accompanying drawing which is a diagrammatic illustration of a medical diagnostic im- 
aging apparatus in accordance with the present invention. 

Referring to the drawing, the apparatus A generates medical diagnostic data which is reconstructed by an 
imager into an electronic image representation. An adaptive filter circuit operates on the electronic image rep- 

10 resentation to improve the image quality and viewabitity thereof. 

Although a magnetic resonance imager is illustrated, the medical diagnostic apparatus A may be a com- 
puterized tomographic scanner, a digital x-ray apparatus, a positron emission scanner, a nuclear camera, or 
other diagnostic apparatus which generates data that is able to be reconstructed into an image representative 
of a region of an examined patient or subject The illustrated magnetic resonance imager includes a field control 

15 means 10 which controls a main, homogeneous polarizing magnetic field through an image region generated 
by electromagnets 12. The field control means 10 also controls gradient magnetic fields created across the 
image region by gradient field coils 14 to provide spatial encoding, phase encoding, and slice select gradients. 
The field control means 10 further generates radio frequency electromagnetic excitation signals which are ap- 
plied to radio frequency coils 16 to excite resonance of dipoles in the image region. A central computer 1 8 con- 

20 trols the relative timing and strengths of the gradient and radio frequency electromagnetic fields. 

Magnetic resonance signals generated by resonating dipoles in the image region are received by the radio 
frequency coils and conveyed to a radio frequency receiver 20. The radio frequency receiver 20 demodulates 
the received signals which are converted from analog to digital format by an A/D converter 22. The digital format 
magnetic resonance medical diagnostic data is supplied to the imager. 

25 The imager under control of the central computer 1 8 reconstructs the medical diagnostic data into an elec- 
tronic image representation. More specifically, the imager reconstructs the diagnostic data which has minor 
noise components into an array of digital pixel values which are degraded by the noise. Each pixel value cor- 
responds to a preselected subregion of the image region, conventionally a corresponding voxel or cubic sub- 
region of an imaged slice. The imager includes a diagnostic data memory 30 which strores the diagnostic data 

30 from the A/D converter 22. The diagnostic data from magnetic resonance echoes are reconstructed into a sin- 
gle image representation. More specific to the preferred embodiment, each echo signal or view is converted 
into a line of ditigal values for a data matrix d(x,y). With Gaussian noise as is found in a magnetic resonance 
imager, each data matrix value is the sum of a data component and a noise component. 

The diagnostic data memory 30 stores the diagnostic data in the data matrix format until a sufficient 

35 amount of data is received for a complex two dimensional Fourier transform means or routine 32 to transform 
the received data into a real image representation for storage in a real image memory 34 and an imaginary 
image representation for storage in an imaginary image memory 36. The sum of a data matrix value and Gaus- 
sian noise, d(x,y)+n g , Fourier transforms linearly into the sum of an image pixel value and Gaussian noise, 
KU) + nV An image select means 38 transfers either real, imaginary, or magnitude pixel values as may be se- 

40 lected by the operator to the filtering means. If the real image is selected, the contents of the real image mem- 
ory 34 are transferred. If the imaginary image is selected, the contents of the imaginary memory 36 are trans- 
ferred. The real and imaginary images again have a Gaussian noise distribution. If a magnitude image is se- 
lected, an array of magnitude pixel values is transferred. Each magnitude pixel value is set equal to the square 
root of the sum of the squares of the corresponding real and imaginary image pixel values, i.e. (I (real) 2 + l(irrv 

45 aginary) 2 )" 2 . The magnitude image has a Rayleigh noise distribution. 

The algorithm implemented by the reconstruction means is, of course, selected in accordance with the 
medical diagnostic apparatus selected. For other diagnostic imagers, other known transform and reconstruc- 
tion techniques are selected. 

The filtering circuit includes an image noise level determining means 40 which calculates a value V(noise) 

so which is indicative of the noise level of the diagnostic data or the overall image. In magnetic resonance and 
many other applications, the diagnostic data stored in memory 30 has substantially signal-independent, un- 
correlated additive noise. This enables image noise statistics to be determined from the data noise statistics. 
In magnetic resonance imaging in which the diagnostic data is a Fourier transform of the image, the noise is 
primarily signal-independent, uncorrected, and additive. 

55 The noise statistics of the data are readily determined from an area of the data which contains no diagnostic 
data or signal, merely noise. In some applications, such a region may be created artificially by appropriate 
placement of the patient, inclusion of phantom or phantom regions, adjustment of the imaging sequence, or 
the like. However, in conventional Fourier transform magnetic resonance imaging, the majority of the signal 



3 



EP0 280 412 B1 



10 



15 



40 



45 



is concentrated in the center of the data matrix d(x,y) stored in memory 30. Accordingly, the probability of find- 
ing a signal-free region is greatest around the periphery of the data matrix. In the preferred embodiment, the 
data around the periphery of the data matrix is examined for uniformity. The most uniform region is assumed 
to be a region with substantially no signal, just noise. The degree of non-uniformity in this region is utilized as 
the indicator of overall image noise. 

A matrix addressing means 42 serially addresses each of a plurality of preselected data matrix points (x,y) 
around the periphery of the data matrix. At each matrix point, a variance means 44 determines the variance 
V(x,y) between the data values corresponding to the addressed matrix point and its neighboring matrix points. 
In the preferred embodiment, the noise at each point V(x,y) is determined from the equation: 



. x+a y+b - 2 



(2), 



where d(f ,k) represents the data values surrounding the matrix point (x,y) and d(x,y) Is the average of data 
values surrounding the data point d(x,y). Each local data variance V(x,y) is compared by a comparing means 
46 with a previously measured variance stored in a memory 48. The smaller of the two variances is returned 

20 to the memory 48 and stored for the next comparison. In this manner, the neighborhood of the data matrix 
with the least signal is determined. For Gaussian noise, the noise variance is equal to the sqaure of standard 
deviation for Gaussian noise oj. 

In the preferred embodiment the data matrix variance determining means 44 includes a means 50 for ser- 
ially receiving the data value from each point surrounding matrix point (x,y), and an averaging means 52 for 

25 determining the average of the neighboring data values. The neighboring points are again addressed serially. 
A subtraction means 54 determines the difference between each neighboring point data value and the average. 
A squaring means 56 squares the difference and a summation means 58 sums the sqaure of the difference 
with the accumulated sum in a cumulative memory 60. The weighting means 62 adjusts the sum by a constant 
to compensate for the size of the neighborhood or the like. 

30 The filtering means further includes a first pixel memory 70 which receives pixel values P(i j) from the re- 
construction means. As discussed above, the pixel values may be the real image pixel values, the imaginary 
image pixel values or the magnitude image pixel values. 

A pixel value average or mean determining means 72 determines the mean or average of pixel values sur- 
rounding each pixel (ij) of the pixel memory 70. In the preferred embodiment, the pixel mean value P(i j) is 

35 determined in accordance with: 

, i+n j+m 



where the values of the constants n and m indicate the size of the neighborhood. 

A pixel value variance determining means 74 determines the variance V(i,j) corresponding to the neigh- 
borhood around each pixel (I j). In the preferred embodiment, the variance is determined in accordance with: 

, 1+n . j4m 2 



l=i-n k=j-ro 



so (4). 



A data noise to image noise converting means 76 converts the Gaussian data noise variance to an image 
noise variance V(noise). For Gaussian image noise, the image noise variance Is equal to the data noise vari- 
ance. For Rayteigh image noise, the image noise variance is equal to the Gaussian data noise variance times 
55 the constant (2-n /2). For other image noise distributions, the data noise may be operated on by other constants 
or functions. 

A weighting factor determining means 78 determines a weighting factor G(i,j) corresponding to each pixel 
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location in accordance with the corresponding pixel variance V(iJ) and the image noise level V(noise). In the 
preferred embodiment, the weighting factor is determined in accordance with: 

G(i.j) = (5). 

v " V(ij) + V(noise) v ' 

5 In most instances, the noise variance V(noise) will be less than or equal to the local variance V(i,j). When 

the noise variance of the image is very low, the weighting factor approaches one. As the noise variance in- 
creases toward the local variance, the weighting factor becomes smaller approaching one half. If the local va- 
riance is less than the noise variance, the weighting factor will become less than one-half and approach zero. 
It must be noted that the noise variance is individually selected for each image and is normally different 

10 from image to image. Thus, the weighting factor at pixels in two images with the same local variance is com- 
monly different due to the different noise variance of the two images. Optionally, other weighting factors which 
weight in accordance with both the local variance and the image noise variance may be selected. 

A mean noise determining means 80 determines the mean magnitude of the image noise correction value 
n. A Gaussian noise distribution has a zero mean magnitude, whereas a Rayleigh distribution does not When 

is a real, imaginary or other image with a Gaussian noise distribution is being processed, the mean noise is zero, 
i.e. n=0. When a magnitude or other image with a Rayleigh noise distribution is being processed, the mean 
noise n is: 

n = (|)i«a 0 

20 Other mean noise values can be calculated for other noise distributions as is known in the art 

Af irst subtraction means 82 subtractively combines each pixel value, P(i,j) with the corresponding average 
neighborhood pixel value, P(i,j). A multiplying means 84 multiplies the difference of the pixel and neighborhood 
values by the weighting factor corresponding to the same pixel G(ij). An adding means 86 adds the correspond- 
ing neighborhood average to the weighted difference. A second subtracting means 88 subtracts the mean 

25 noise correction value from the sum. This sum is stored in a filtered image memory 90 and displayed on a 
video monitor or other display means 92. That is, each filtered pixel value P*(i,j) is equal to: 

P'ftP* POD* ^"^-n (7). 
\w i-vdJ V0j) + V(noise) v ' 

With this relationship, with less data noise variance, the corresponding pixel value P(i,j) is weighted most 

30 heavily. Similarly, in a local region with alot of variance, the actual pixel value is again weighted heavily so as 

not to blur the iocai resolution. However, in images that are very noisy or in regions that are very uniform, the 

neighborhood average is weighted more heavily. 

The invention has been described with reference to the preferred embodiment. Obviously, modifications 

and alterations will occur to others upon reading and understanding the preceding detailed description. It is 

35 intended that the invention be construed as including all such alterations and modifications insofar as they 

come within the scope of the appended claims or the equivalents thereof. 



Claims 

40 

1. A method of imaging comprising: collecting image data; converting the collected data into an electronic 
image representation which includes an array of pixel values; for each pixel value determining (i) a vari- 
ance between the pixel value and neighboring pixel values and (ii) an average of the neighboring pixel 
values; and replacing each pixel value of the array by a weighted combination of the said pixel value and 

45 the average of pixel values neighboring said pixel value; characterised in that the method further includes 

determining noise statistics of the collected data; and in that the weighting is in accordance with the de- 
termined noise statistics and the variance between said pixel value and its neighboring pixel values. 

2. A method according to Claim 1 wherein the noise statistics of the collected data are determined by: digi- 
50 tizing the collected data and storing the digitized data in a data memory; determining a data variance be- 
tween each of a plurality of data values stored in the data memory and its neighboring data values; and 
comparing each determined data variance with previously determined data variances to determine a mini- 
mum data variance, the minimum data variance being indicative of the noise statistics, whereby the noise 
variance of the collected diagnostic data is determined. 

55 3. A method according to Claim 2 wherein the data variance is determined only for selected data values 
whose surrounding data values have a higher probability of representing substantially no signal such that 
the data variance is primarily indicative of noise variance rather than imaged structure. 
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4. A method according to Claim 2 or Claim 3 further including: generating the data with a magnetic resonance 
imager; organizing the digitized data values in the data memory in a rectangular matrix array; and, in the 
noise variance determination, determining the variance of data values disposed generally around the per- 
iphery of the matrix. 

5. A method according to any one of the preceding claims wherein the weighting is proportional to the ratio 
of the corresponding data value variance to a sum of the corresponding pixel value variance and the de- 
termined noise statistics. 

6. A method according to any one of the preceding claims wherein the replacing step includes subtractively 
combining replaced pixel values with the average of neighboring pixel values and multiplying the differ- 
ence with a weighting function, the weighting function being determined in accordance with the deter- 
mined noise statistics and the pixel value variance between the replaced pixel and its neighboring pixel 
values. 

7. A method according to Claim 6 wherein the weighting function is proportional to a ratio of the pixel value 
variance corresponding to the replaced pixel and the determined noise statistics. 

8. A method according to Claim 7 wherein the replacing step further includes summing the average of pixel 
values neighboring the replaced pixel value with the product of said difference and said weighting func- 
tion. 

9. A method according to Claim 8 further including: determining a mean noise level value from the deter- 
mined noise statistics; and wherein the replacing step further includes subtracting the mean noise level 
value from said weighted combination. 

10. An imaging apparatus comprising: data generating apparatus (A) for generating image data indicative of 
at least a selected region of a body; a data memory means (30) for storing image data from the data gen- 
erating apparatus (A); a transform means (32, 34, 36) for converting the data from the data memory 
means (30) into an electronic image representation which includes an array of pixel values; an average 
pixel value determining means (38, 70, 72) for determining an average pixel value of pixel values neigh- 
boring each pixel value; and combining means (82, 84, 86, 88) for adjusting each pixel value, in accordance 
with the corresponding said average pixel value; characterised in that the apparatus further includes a 
data noise statistic determining means (40, 76) for determining a noise statistic from the data stored in 
said data memory means (30); and in that said combining means (82, 84, 86, 88) combines each pixel 
value with the corresponding average pixel value and the data noise statistic. 

11. An apparatus according to Claim 10 further including: 

a variance determining means (74) for determining a variance between each pixel value and its neigh- 
boring pixel values; a weighting function determining means (78) for determining a weighting function for 
each pixel value in accordance with (i) a data noise level from the noise statistic determining means (40, 
76) and (ii) the variance corresponding to the same pixel value; and wherein the combining means (82, 
84, 86, 68} weights a combination of each pixel value and the corresponding average pixel value with the 
weighting function. 

12. An apparatus according to Claim 10 or Claim 12 wherein the combining means (82, 84, 86, 88) includes: 
a subtraction means (82) for subtractively combining each pixel value with its corresponding average of 
neighboring pixel values; means for weighting (84) the difference from the subtractive combining means 
with the weighting factor; and an adding means (86) for combining the corresponding average of neigh- 
boring pixel values with the weighted difference. 

13. An apparatus according to any one of Claims 10 to 12 further including: an analog-to-digital converter 
(22) for digitizing the image data; and wherein the data memory means (30) stores the digitized data in 
a rectangular matrix array. 

14. An apparatus according to Claim 13 wherein the image data generating apparatus (A) Is a magnetic res- 
onance imager (A); wherein the digital data values around a periphery of the rectangular matrix array tend 
to have a smaller signal contribution; and further including: a data memory addressing means (42) for 
addressing the data values around the periphery. of the matrix; a data variance determining means (44) 
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for determining a data value variance between each addressed matrix value and its neighboring matrix 
values; and a comparing means (46, 48) for comparing each determined data value variance of previously 
addressed data points, the smallest of the determined data value variances being the data noise statistic, 
the comparing means (46, 48) being operatively connected with' the combining means (82, 84, 86, 88) for 
5 supplying the smallest data value variance thereto. 

1 5. An apparatus according to any one of Claims 10 to 14 further including: a mean image noise determining 
means (80) for determining a mean image noise value from the output of the data noise statistic means 
(40, 76); and wherein the combining means (82, 84, 86, 88) subtracts the mean image noise value from 

to the combination of each pixel value and the corresponding average pixel value. 

1 6. An apparatus according to any one of Claims 10 to 15 further including an image memory (90) connected 
with the combining means (82, 84, 86, 88) for storing the adjusted pixel values produced thereby. 



Paten tansp ruche 

1. Abbildungsverfahren, umfassend: Erfassen von Bilddaten; Umsetzen der erfa&ten Daten in eine elek- 
tronische Bi Id re presentation, die ein Feld aus Pixelwerten umfa&t; Bestimmen fur jeden Pixelwert eine 

20 Varianz (i) zwischen dem Pixelwert und benachbarten Pixelwerten und (ii) einen Durchschnittswert der 

benachbarten Pixelwerte; und Ersetzen jedes Pixelwerts des Feldes durch eine gewichtete Kombination 
dieses Pixelwerts und des Durchschnitts von Pixelwerten, die diesem Pixelwert benachbart sind; 
dadurch gekennzeichnet, 

dad das Verfahren ferner die Ermtttlung von Rauschstatistiken der erfa&ten Daten umfa&t; und daft die 
25 Wichtung gemafc den ermittelten Rauschstatistiken und der Varianz zwischen dem Pixelwert und dessen 

Nachbarpixelwerten ist. 

2. Verfahren nach Anspruch 1, in welchem die Rauschstatistiken der erfa&ten Daten bestimmt werden 
durch: Digitalisieren der erfaftten Daten und Speichern der digital isierten Daten in einem Datenspeicher; 
Ermitteln einer Datenvarianz zwischen jedem mehrerer Datenwerte, die im Datenspeicher gespeichert 
sind, und dessen benachbarten Datenwerten; und Vergieichen jeder ermitteiten Datenvarianz mitzuvor 
ermittelten Datenvarianzen zur Ermittlung einer minimalen Datenvarianz, wobei die minimale Datenva- 
rianz die Rauschstatistiken anzeigt, wodurch die Rauschvarianz der erfa&ten Diagnosedaten ermittelt 
wird. 

35 3. Verfahren nach Anspruch 2, in welchem die Datenvarianz nur fur selektierte Datenwerte ermittelt wird, 
deren umgebende Datenwerte eine hohere Wahrscheinlichkeit aufweisen, im wesentlichen kein Signal 
darzustellen, derart, da& die Datenvarianz primar eine Rauschvarianz statt eine abgebildete Struktur an- 
zeigt 

40 4. Verfahren nach Anspruch 2 oder 3, ferner umfassend: Erzeugen der Daten mit einem Magnetresonanz- 
Abbildungsgerat; Vornehmen eines organisatorischen Aufbaus der digitafisierten Datenwerte im Daten- 
speicher in einer rechteckigen Matrixanordnung; und Ermitteln der Varianz von Datenwerten bei der 
Rauschvarianzermittlung, die im wesentlichen urn den Umfangsbereich der Matrix angeordnet sind. 

45 5. Verfahren nach einem der vorhergehenden AnsprOche, in welchem die Wichtung proportional zum Ver- 
haltnis der entsprechenden Datenwertvarianz zu einer Summe derentsprechenden Pixelwertvarlanz und 
der ermittelten Rauschstatistiken ist 

6. Verfahren nach einem der vorhergehenden AnsprQche, in welchem der ersetzende Schritt die subtraktive 
so Kombination ersetzter Pixelwerte mit dem Durchschnittswert benachbarter Pixelwerte und die Multipli- 

kation der Differenz mit einer Wichtungsf unktion umfa&t, wobei die Wichtungsf unktion entsprechend den 
ermittelten Rauschstatistiken und der Pixelwertvarianz zwischen dem ersetzten Pixel und dessen Nach- 
barpixelwerten ermittelt wird. 

65 7. Verfahren nach Anspruch 6, in welchem die Wichtungsfunktion proportional zu einem Verhiltnis der Pi- 
xelwertvarianz entsprechend dem ersetzten Pixel und den ermittelten Rauschstatistiken ist 

8. Verfahren nach Anspruch 7, in welchem der ersetzende Schrittferner die Summierung des Durchschnitts- 
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werts von Pixelwerten, die dem ersetzten Pixelwert benachbart sind, zum Produkt der Differenz und der 
Wichtungsfunktion umfaBL 

9. Verfahren nach Anspruch 8 ( ferner aufweisend: Ermitteln eines mittleren Rauschpegelwerts aus den er- 
5 mittelten Rauschstatistiken; und in welcbem der ersetzende Schritt ferner die Subtraktion des mittleren 

Rauschpegelwerts von der gewichteten Kombination umfa&t 

1 0. Abbildungsvorrichtung, aufweisend: eine Datenerzeugungsvorrichtung (A) zum Erzeugen von Bilddaten, 
die zumindest eine selektierte Region eines Kdrpers anzeigen; eine Datenspeichereinrichtung (30) zum 

10 Speichern von Bilddaten von der Datenerzeugungsvorrichtung (A); eine Transformationseinrichtung (32, 

34, 36) zum Umsetzen der Daten aus der Datenspeichereinrichtung (30) in eine elektronische Bildre pre- 
sentation, die ein Feld aus Pixetwerten umfa&t; eine Durchschnittspixeiwert-Ermittlungseinrichtung (38, 
70, 72) zum Ermitteln eines durchschnittlichen Pixelwerts von Pixelwerten, die jedem Pixelwert benach- 
bart sind; und eine Kombinationseinrichtung (82, 84, 86, 88) zur Einstetlung Jedes Pixelwerts gemaa dem 

15 entsprechenden solchen Durchschnittspixelwert; 

dadurch gekennzelchnet, 

daft die Vorrichtung ferner eine Datenrauschstatistik-Ermittlungseinrichtung (40, 76) aufweist, die eine 
Rauschstatlsti k aus den in der Speichereinrichtung (30) gespeicherten Daten ermittelt; und daS die Kom- 
binationseinrichtung (82, 84,86, 88) jeden Pixelwert m it dem entsprechenden Durchschnittspixelwert und 
OA der Datenrauschstatisti k kombiniert. 

20 

11. Vorrichtung nach Anspruch 10, ferner aufweisend: eine Varianzermittlungseinrichtung (74) zum Ermitteln 
einer Varianz zwischen jedem Pixelwert und dessen benachbarten Pixelwerten; eine Wichtungsfunktions- 
Ermittlungseinrichtung (78) zum Ermitteln einer Wichtungsfunktion fur jeden Pixelwert gemafc (i) einem 
Datenrauschpegel von der Rauschstatistik-Ermittlungseinrichtung (40, 76) und (ii) der Varianz entspre- 
chend demselben Pixelwert; und in welcher die Kombinationseinrichtung (82, 84, 86, 88) eine Kombina- 
tion jedes Pixelwerts und des entsprechenden Durchschnittspixelwerts mit der Wichtungsfunktion 
wichtet. 

1 2. Vorrichtung nach Anspruch 1 0 oder Anspruch 1 2, in welcher die Kombinationseinrichtung (82, 84, 86, 88) 
30 umfaat eine Subtraktionseinrichtung (82) zur subtraktiven Kombination jedes Pixelwerts mit seinem ent- 
sprechenden Durchschnittswert von benachbarten Pixelwerten; eine Einrichtung zum Wichten (84) der 
Differenz von der subtraktiven Kombinationseinrichtung mit dem Wichtungsfaktor; und eine Additions- 
einrichtung (86) zur Kombination des entsprechenden Durchschnittswerts benachbarter Pixelwerte mit 
der gewichteten Differenz. 

35 

13. Vorrichtung nach einem der Ansprflche 10 bis 12, ferner aufweisend: einen Analog/Digital-Wand ler (22) 
zur Digitalisierung der Bilddaten; und in welcher die Datenspeichereinrichtung (30) die digitalisierten Da- 
ten in einem rechteckigen Matrixfeld speichert. 

ao 14. Vorrichtung nach Anspruch 13, in welcher die Bilddatenerzeugungseinrichtung (A) ein Magnetresonanz- 
Abbildungsgerat (A) 1st; in welcher die digitalen Datenwerte urn einen Umfang des rechteckigen Matrix- 
feldes dazu neigen, einen geringeren Signal beitrag aufzuweisen; und ferner aufweisend: eine 
Datenspeicheradressierungseinrichtung (42) zum Adressieren der Datenwerte urn den Umfang der Ma- 
trix; eine Datenvarianzermittlungseinrichtung (44) zum Ermitteln einer Datenwertvarianz zwischen jedem 

45 adressierten Matrixwert und seinen benachbarten Matrixwerten; und eine Vergleichseinrichtung (46, 48) 

zum Vergleichen jeder ermittelten Datenwertvarianz zuvor adressierter Datenpunkte, wobei die geringste 
der ermittelten Datenwertvarianzen die Datenrauschstatisti k (KenngroBe) ist, wobei die Vergleichsein- 
richtung (46, 48) betrlebswirksam mit der Kombinationseinrichtung (82, 84, 86, 88), um dieser die gering- 
ste Datenwertvarianz zuzuf uhren, verbunden ist. 

so 

15. Vorrichtung nach etnem der Ansprflche 10 bis 14, ferner aufweisend: eine Bildrauschmittelwert-Ermitt- 
lungseinrichtung (80) zum Ermitteln eines mittleren Bildrauschwerts aus dem Ausgangssignal der Da- 
tenrauschstatistikeinrichtung (40, 76); und in welcher die Kombinationseinrichtung (82, 84, 86, 88) den 
mittleren Bildrauschwert von der Kombination jedes Pixelwerts und des entsprechenden 

55 Durchsch nits pixelwerts subtrahlert. 

16. Vorrichtung nach einem der Ansprflche 10 bis 15, ferner aufweisend einen Bildspeicher (90), der mit der 
Kombinationseinrichtung (82, 84, 86, 88) verbunden ist und die hierdurch erzeugten eingestellten Pixel- 
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werte speichert. 



Revendlcations 

5 

1. Procede de formation d'image, comprenant la collecte de donnees d' image, ta conversion des donnees 
d' image en une representation electron ique d'image qui comprend une matrice de valeurs d'elements 
d'image, et, pour chaque valeur d'element d'image, la determination (i) d'une variance entre la valeur de 
I'element d'image et les valeurs des elements d'image voisins, et (ii) d'une moyenne des valeurs des ele- 

10 ments d'image voisins, et le remplacement de chaque valeur d'element d'image de la matrice par une 

comblnaison ponderee de la valeur de I'element d'image et de la moyenne des valeurs d'elements d'image 
voisins de la valeur de ('element d'image, caracterise en ce que le precede comporte en outre la deter- 
mination de la statistique de bruit des donnees collectees, et en ce que la ponderation s'effectue sutvant 
la statistique de bruit determinee et la variance entre la valeur de I'element d'image et les valeurs des 

15 elements d'image voisins de celuf-ci. 

2. Precede selon la revendication 1, dans lequel la statistique de bruit des donnees collectees est determi- 
nee par numerisation des donnees collectees et memorisation des donnees numerisees dans une me- 
moire de donnees, determination d'une variance des donnees entre chacune de plusieurs valeurs de don- 

20 nees conservees dans la memoire de donnees et les valeurs des donnees voisines, et comparaison de 

chaque variance des donnees determinee a des variances de donnees determines precedemment pour 
la determination d'une variance minimale des donnees, la variance minimaie des donnees etant repre- 
sentative de la statistique de bruit, si bien que la variance du bruit des donnees collectees de diagnostic 
est determinee. 

25 

3. Procede selon la revendication 2, dans lequel la variance des donnees n'est determinee que pour des 
valeurs choisies de donnees dont les valeurs des donnees environnantes ont une grande probabilite de 
representer pratiquement I'absence de signal, si bien que la variance des donnees est essentiellement 
representative de la variance du bruit plutdt que de ia structure de I' image. 

30 

4. Procede selon la revendication 2 ou 3, comprenant en outre la creation de donnees avec un appareil de 
formation d'image par resonance magnetique, I 'organisation des valeurs numerisees de donnees dans 
la memoire de donnees sous forme d'une matrice rectangulaire et, au cours de la determination de la va- 
riance de bruit, la determination de la variance des valeurs de donnees disposees de facon generate a 

^ la peripheric de la matrice. 

5. Procede selon Tune quelconque des revendications precedentes, dans lequel la ponderation est propor- 
tion n el le au rapport de la variance de la valeur des donnees correspond antes et d'une somme de la va- 
riance de la valeur de ('element d'image correspondant et de la statistique de bruit qui est determinee. 

40 6. Procede selon Tune quelconque des revendications precedentes, dans lequel I'etape de remplacement 
comprend ta combinaison par soustraction des valeurs des elements d'image qui sont remplacees par 
la moyenne des valeurs d'elements d'image voisins et par multiplication de la difference par une fonction 
de ponderation, la fonction de ponderation etant determinee d'apres la statistique de bruit determinee et 
la variance de la valeur de I'element d'image entre I'element d'image remplace et les valeurs des elements 

45 d'image voisins. 

7. Procede selon la revendication 6, dans lequel la fonction de ponderation est proportionnetle a un rapport 
de la variance de la valeur de I'element d'image correspondant a I'element d'image remplace et de la sta- 
tistique de bruit determinee. 

50 

8. Procede selon la revendication 7, dans lequel I'etape de remplacement comporte en outre la sommation 
de la moyenne des valeurs d'elements d'image voisins de la valeur de I'element d'image remplace par 
le produit de la difference et de la fonction de ponderation. 

55 9. Procede selon la revendication 6, comprenant en outre la determination de la valeur d'un niveau moyen 
de bruit tiree de la statistique de bruit determin6e, et dans lequel I'etape de remplacement comporte en 
outre la soustraction de la valeur du niveau moyen de bruit de la combinaison ponderee. 
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10. Appareil de formation d'image comprenant un appareil (A) generateur de donnees d' image representa- 
tives d'au moins une region choisie d'un corps, un dis posit if (30) a memoire de donnees destine a conser- 
verles donnees d'image provenantde I'appareil (A) generateur de donnees, un disposittf (32, 34, 36) de 
transformation qui convertit les donnees du dispositif (30) a memoire de donnees en une representation 

s electronique d'image qui comprend une mat rice de valeurs d'elements d'image, un dispositif (38, 70, 72) 

de determination d'une valeur moyenne d'element d'image des valeurs des elements d'image voislns de 
chaque valeur d'element d'image, et un dispositif (82, 84, 86, 88) de combinaison destine a ajuster chaque 
valeur d'element d'image d'apres la valeur moyenne correspondante d'element d'image, caracterfse en 
ce que I'appareil comporte en outre un dispositif (40, 76) de determination d'une statistique de bruit des 

10 donnees destine a determiner une statistique de bruit a partir des donnees conservees dans le dispositif 

(30) a memoire de donnees, et en ce que le dispositif de combinaison (82, 84, 86, 88) combine chaque 
valeur d'element d'image a la valeur moyenne correspondante d'element d'image et a la statistique de 
bruit des donnees. 



15 11. Appareil selon la revendication 10, comprenant en outre un dispositif (74) de determination d'une variance 
entre chaque valeur d'element d'image et les valeurs des elements d'image voisins, un dispositif (78) de 
determination d'une fonction de ponderation pour chaque valeur d'element d'image d'apres (i) un niveau 
de bruit de donnees provenant du dispositif (40, 76) de determination d'une statistique de bruit, et (ii) la 
variance correspondent a la meme valeur d'element d'image, et dans lequel le dispositif de combinaison 

20 (82, 84, 86, 88) pondere une combinaison de chaque valeur d'element d'image et de la valeur moyenne 

correspondante d'element d'image a la fonction de ponderation. 



12. Appareil selon la revendication 10 ou 11, dans lequel le dispositif de combinaison (82, 84, 86, 88) 
comprend un dispositif (82) de soustraction destine a combiner de maniere soustractive chaque valeur 
d'element d'image et la moyenne correspondante des valeurs des elements d'image voisins, un dispositif 
(84) de ponderation de la difference provenant du dispositif de combinaison par soustraction a I'aide du 
facteur de ponderation, et un dispositif addittonneur (86) destine a combiner la moyenne correspondante 
des valeurs des elements d'image voisins a la difference ponderee. 

13. Appareil selon Tune quelconque des revendications 1 0 a 12, comprenant en outre un convertisseur ana- 
log ique-numerique (22) destine a numeriser les donnees d'image, et dans lequel le dispositif (30) a me- 
moire de donnees conserve les donnees numerisees sous forme d'une matrice rectangulaire. 



14. Appareil selon la revendication 1 3, dans lequel I'appareil (A) generateur des donnees d'image est un ap- 
pareil (A) de formation d'image par resonance magnetique, dans lequel les valeurs des donnees nume- 

35 riques placees a la peripheric de la matrice rectangulaire ont tendance a presenter une plus petite contri- 

bution au signal, et comprenant en outre un dispositif (42) d'adressage de memoire de donnees destine 
a adresser les valeurs des donnees a la peripheric de la matrice, un dispositif (44) de determination d'une 
variance de la valeur des donnees entre chaque valeur adressee de la matrice et les valeurs voisines de 
la matrice, et un dispositif (46, 48) de comparaison de chaque variance determinee pour une valeur de 

40 donnees pour les points de donnees anterieurement adresses, la plus petite des variances determinee 

pour les valeurs de donnees etant la statistique de bruit des donnees, le dispositif de comparaison (46, 
48) etant connecte pendant le fonctionnement au dispositif de combinaison (62, 84, 86, 88) af in qu'il trans- 
mette la plus petite variance de valeur de donnees au dispositif de comparaison. 

45 15. Appareil selon I'une quelconque des revendications 10 a 14, comprenant en outre un dispositif (80) de 
determination de bruit moyen de I'image destine a determiner une valeur moyenne du bruit de I'image a 
partir du signal de sortie du dispositif (40, 76) de determination de statistique de bruit des donnees, et 
dans lequel le dispositif de combinaison (82, 84, 86, 88) soustrait la valeur moyenne du bruit de I'image 
de la combinaison de chaque valeur d'element d'image et de la valeur moyenne correspondante d'element 

50 d'image. 

16. Appareil selon I'une quelconque des revendications 10 a 15, comprenant en outre une memoire (90) 
d'image connectee au dispositif de combinaison (82, 84, 86, 88) et destinee a memoriser les valeurs ajus- 
tees d 'elements d'image ainsi produites. 

55 
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